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Abstract: We present "Before and After" global oscillation solutions, as well as 
predicted "Before and After" values and ranges for ten future solar neutrino observ- 
ables (for BOREXINO, KamLAND, SNO, and a generic p — p neutrino detector). 
The "Before" case includes all solar neutrino data (and some theoretical improve- 
ments) available prior to April 20, 2002 and the "After" case includes, in addition, 
the new SNO data on the CC, NC, and day-night asymmetry. We have performed 
global analyses using the full SNO day-night energy spectrum and, alternatively, 
using just the SNO NC and CC rates and the day-night asymmetry. The LMA solu- 
tion is the only currently allowed MSW oscillation solution at ~ 99% CL. The LOW 
solution is allowed only at more than 2.5a, SMA is now excluded at 3.7a or 4.7cr 
depending upon analysis strategy, and pure sterile oscillations are excluded at more 
than 4.7cr. Small mixing angles are "out" (pure sterile is "way out"); MSW with 
large mixing angles is definitely "in." Vacuum oscillations are allowed at 3a, but not 
at 2a. Precise maximal mixing is excluded at 3.2a for MSW solutions and at more 
than 2.8a for vacuum solutions. Most of the predicted values for future observables 
for the BOREXINO, KamLAND, and future SNO measurements are changed only 
by minor amounts by the inclusion of the recent SNO data. In order to test the 
robustness of the allowed neutrino oscillation regions that are inferred from the mea- 
surements and the predicted values for future solar neutrino observables, we have 
carried out calculations using a variety of strategies for analyzing the SNO and other 
experimental data. 
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1. Introduction 

The goal of this paper is to assess the impact of recent SNO measurements % § 
on the allowed regions of neutrino oscillation parameters and on the predicted values 
of the most important future solar neutrino observables. The SNO collaboration 
has reported a neutral current (NC) measurement of the active 8 B solar neutrino 
flux and related measurements of the day-night asymmetry, as well as improved 
determinations of the charged current (CC) and neutrino-electron scattering rate. 
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We are concerned that global solutions for oscillation parameters depend upon 
the assumption that the errors are well understood for all of the reported measure- 
ments in the chlorine, gallium, Super-Kamiokande, and SNO experiments. There 
are many cases in the past for which similar assumptions have proved misleading. 
Therefore, we focus our study on determining the robustness of our conclusions re- 
garding the currently allowed oscillation parameters and the predicted values of new 
solar neutrino observables. We test the robustness of the conclusions about allowed 
oscillation parameters by using three different analysis strategies. We also compute 
predicted values for future solar neutrino observables with the oscillation parameters 
that are currently allowed as well with the oscillation parameters that were allowed 
before the recent SNO measurements. In addition, we treat the SNO data in two 
different ways: a) with the aid of a two-step process (discussed in section |2~5|) and 
b) using the full SNO day-night energy spectrum (see section |5.4|) used by the SNO 
collaboration. In our view, a necessary condition for a result to be regarded as "ro- 
bust" is that the result not change significantly as we vary the analysis procedures 
among the different plausible possibilities listed above. 

Perhaps the most remarkable and encouraging result of our analysis is that the 
allowed regions for the oscillation parameters, Am 2 and tan 2 9, the predicted values 
of ten future solar neutrino observables, and the inferred total 8 B neutrino flux are 
all rather robust with respect to the choices among the different analysis procedures. 
The inferences from the available data are relatively independent of the details of 
the analysis procedures. This conclusion will be justified quantitatively in section |5] 
and summarized in section |6]. 

We begin in section |2| by deriving, using the two-step procedure for the SNO data, 
the currently allowed regions in neutrino oscillation space that are obtained with 
three different analysis strategies (see figure P, each strategy previously advocated 
by a different set of authors. The global solutions obtained here are calculated 
using the methods described in our recent paper [Q (see especially section 3.3 of 
ref. ||) but also include some refinements in addition to the new SNO data. For 
example, we take account of the energy dependence and correlations of the errors 
in the neutrino absorption cross sections for the chlorine and gallium solar neutrino 
experiments as described in the Appendix of ref. ]5j]. We also include the recently 
reported SAGE || data for 11 years of observation and the zenith angle-recoil energy 
spectrum data presented by the Super-Kamiokande collaboration after 1496 days of 
observations Q] . Where required, we use the predicted fluxes and their errors from 
the BPOO standard solar model ||. 

We present in section |3], and especially in figure [3], a "Before and After" com- 
parison of the globally allowed neutrino oscillation solutions (see figure ^|). In the 
"Before" case, we use all the solar neutrino data (see refs. || |K], [11], [12], [13], [L4], [T^j ) 
that were published or had appeared publicly before April 20, 2002, the date that 
the SNO NC and day-night asymmetry were first published. In the "After" case, 
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we include in addition measurements reported in the two recent papers [|T], |] by the 
Sudbury Neutrino Observatory (SNO) collaboration. 

What are the predicted values of the ten most informative quantities that can 
be measured in the reasonably near future in solar neutrino experiments? We use in 
section |] the allowed regions in neutrino parameter space (obtained in section |2|) to 
predict in table |2| the expected range of the most promising quantities that can be 
measured accurately in the BOREXINO |TE[ and KamLAND |T7[ 7 Be solar neutrino 



experiments and in the KamLAND reactor experiment, as well as the spectrum 
distortion and the day-night asymmetry in the SNO CC measurements. We also 
include predictions for a generic p — p neutrino-electron scattering detector. To 
assess the robustness of the predictions, we compare the values predicted using the 
"After April 20, 2002" global solution (table ||) with the values predicted using the 
"Before April 20, 2002" solution (table g). 

We summarize and discuss our conclusions in section 

We do not discuss in this paper the implications of the agreement between the 
measured [1|] flux of active 8 B solar neutrinos and the predicted || standard solar 
model 8 B neutrino flux. The agreement is accidentally too good to be true [see 
eq. ( p.6|) ]. As more measurements are made of the neutrino flux and of the solar 
model parameters the agreement should become less precise. We are aware of several 
recent and ongoing measurements, which are currently not in good agreement, for the 
low energy cross section factor Sn, to which the calculated standard solar model 8 B 
neutrino flux is proportional. Until the new laboratory measurements of Si 7 converge 
to a better defined range, we continue to use the standard value adopted in BP00 2 . 

2. Global oscillation solutions 



We describe in section [2.1| the oscillation solutions that are allowed with three dif- 
ferent analysis prescriptions (see figure [1]). We have at different times used all three 

1 Several papers [u^, [2(| have appeared essentially contemporaneously with the present paper 
and treat some of the same topics with somewhat similar results, although refs. jl8]. have not 

calculated predictions for the ten future solar neutrino observables studied in the present paper. On 
a technical level, as far as we can tell, these papers have not included the correlations and the energy 
dependences of the neutrino absorption cross sections for the chlorine and gallium experiments (see 
the Appendix of ref. H). Also, the potential effects of distortions on the interpretation of the SNO 



data in terms of individual rates and their error correlations (see section 2J3) were not treated in 
detail in the originally-posted versions, although more complete treatments have been made in later 
versions [Q. Both of these effects are included in the present paper. A concise but insightful and 
informative discussion of the effects of the recent SNO measurements on solar neutrino oscillations 
is given in the original SNO NC paper H. The interested reader may wish to consult in addition 
a number of recent papers, refs. p2|T [23|, gjj , p5| p6| , that have determined from a variety of 
perspectives the allowed solar neutrino oscillation solutions following the June, 2001 announcement 
of the SNO CC measurement ||. 

2 For an insightful discussion of the predicted and measured total 8 B neutrino flux, see ref. |18|. 
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of the analysis strategies and various colleagues have advocated strongly one or the 
other of the strategies described here. Given the recent SNO NC measurement, we 
now prefer the strategy in which the 8 B neutrino flux is treated as a free parameter. 
This strategy is implemented in figure 0a (see the discussion below). A comparison 
of the results obtained using the three strategies allows one to test the robustness of 
any conclusion to the method of analysis. 

We present in section 272 and table [T] the best-fit oscillation parameters for 
the allowed and the disfavored solutions, treating the 8 B neutrino flux as a free 
parameter. We discuss in this section the CL at which different oscillation solutions 
are acceptable. 

In section |2.3| , we present and discuss the allowed ranges for Am 2 , tan 2 9, and 
the total active flux of 8 B solar neutrinos. 



We describe in section |2]<| and in figure the predicted dependence of the survival 
probability as a function of energy and of day or night for the best-fit LMA, LOW, 
and vacuum solutions. Figure ^ in particular provides a succinct overview of the 
energy and day-night dependences of the best-fit survival probabilities. 

The SNO experiment detects CC, NC, ES [y — e scattering), and background 
events. The rates from these different processes are correlated because they are, 
with the present data, observed most accurately in a mode in which all of the events 
are considered together and a simultaneous solution is made for each of the separate 
processes using their known angular dependences (with respect to the solar direction) 
and their radial dependences in the detector, as well as information from direct 
measurements of the background [|T], |2], U . Details of how this analysis is done are 
given in refs. ^|, |3|. 

Since the different processes are coupled together in the analysis, the inferred 
values for the measured fluxes in each of the CC, NC, and ES modes depend upon 
the assumed distortion of the CC and ES recoil energy spectra, which in turn depend 
upon the assumed Am 2 and tan 2 9. In order to avoid this cycle, the SNO collab- 
oration presented results for the CC, NC, and ES fluxes that were determined 
by assuming that the CC and ES recoil energy spectra are undistorted by neutrino 
oscillations or any other new physics. 

The fluxes inferred using the hypothesis of undistorted energy spectra has been 
used by ref. |2(| in their analyses of the SNO data. This approximation is excellent 
for the LMA and LOW solutions, as we shall show later in this section (see, e.g., 
results quoted in section ^2|) , because for these solutions the expected distortions 
are indeed small. The approximation is less accurate for other solutions in which the 
distortions are more significant. 

In order to obtain the results presented in the following subsections, we used 
a two-step procedure to take account, for each value of Am 2 and tan 2 9, of the 
energy dependence of the survival probability. We have used in these calculations 
experimental data provided by the SNO collaboration B. We describe how we have 
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Figure 1: Global neutrino oscillation solutions for three different analysis 
strategies. The strategy used in constructing panel (a) treats the 8 B solar neutrino 
flux as a free parameter to be determined by the experimental data together with Am 2 
and tan 2 9. The strategies corresponding to panels (b) and (c) include the theoretical un- 
certainty in the 8 B neutrino flux, but differ in how they treat the total rate measured in 
the Super-Kamiokande experiment. The input data used in constructing figure [l] include 
the neutrino fluxes and uncertainties predicted by the BPOO solar model |8| and the to- 
tal measured CC and NC event rates from the SNO experiment [Q], the SNO day-night 
asymmetry Q, the Chlorine [ 10 1 and Gallium (average d) |, 0, U, HI event rates, as well 



as the zenith angle-recoil energy spectrum data presented by Super-Kamiokande 14]. The 
rates from the GALLEX/GNO and SAGE experiments have been averaged to provide a 
unique data point (72.4 ± 4.7 SNU). The CL contours shown in the figure are 90%, 95%, 
Vo, and 99.73% (3a). The global best-fit points are marked by a star. 



carried out the details of this analysis in section p75| . The results obtained using the 
full SNO day-night energy spectrum are presented in section |[ 

2.1 Three strategies 

Figure [l| shows the allowed ranges of the neutrino oscillation parameters, Am 2 and 
tan 2 6, that were computed using the three different analysis approaches that have 
been used previously in the literature. We use the analysis methods and procedures 
described in refs. [|, [|, [27], |29|, 0), see especially section 3.3 of ref. || and 



the Appendix of ref. || . We follow refs. [51], |32| in using tan 2 9 (rather than sin 2 26) 



in order to display conveniently the solutions on both sides of 9 = 7r/4. 

Figure |l|a presents the result for our standard analysis, (a), the free 8 B analysis. 
The strategy used in this standard analysis takes account of the BPOO predicted fluxes 
and uncertainties for all neutrino sources except for 8 B neutrinos. The normalization 
of the 8 B neutrino flux is treated as a free parameter for analysis strategy (a) but 
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is constrained by the BPOO prediction and error for strategies (b) and (c). Analysis 
strategy (a) considers all the experimental data except for the Super-Kamiokande 
total event rate. The recoil electron zenith angle-recoil energy spectrum data repre- 
sent the Super-Kamiokande total rate in this approach. Details of the treatment of 
the Super-Kamiokande errors are given in Sec. |5.4| . 

Figure |I|b displays the results of a calculation, analysis strategy (b), which is 
the same as for the standard case, figure |T]a, except that the 8 B neutrino flux is 
constrained by the BPOO solar model prediction. Prior to the SNO measurement 
of the NC flux fl[], (b) was our standard analysis strategy (cf. ref. ||). Following 
the SNO NC measurement, we prefer to use the neutrino data to determine the flux 
normalization and therefore to test, rather than assume, the standard solar model 
prediction for the 8 B neutrino flux. 

Figure |TJc was constructed by an analysis similar to that used to construct fig- 
ure |I|b except that for figure 0b the total Super-Kamiokande rate is included explicitly 
together with a free normalization factor for the zenith angle-recoil energy spectrum 
of the recoil electrons. This procedure has been used especially effectively by the 
Bari group [^3 



In all cases, we have accounted for the energy and time dependence of the Super- 
Kamiokande data by using their zenith angle-recoil energy spectrum data. The 
available Super-Kamiokande zenith angle-recoil energy spectrum consists of 44 data 
points, corresponding to six night bins and one day bin for six energy bins between 
5.5 and 16 MeV electron recoil energy, plus two daily averaged points for the lowest 
(5.0 < E < 5.5 MeV) and the highest (E > 16 MeV) energy bins. Alternatively, 
one could use their day-night energy spectra as given in 19 energy bins each for the 
day and for the night periods. Using the more complete zenith angle-recoil energy 
spectrum data allows for a better discrimination between the LMA and the LOW 
solutions. Within the LMA regime, oscillations in the Earth are rapid and therefore 
LMA predicts a rather flat distribution in zenith angle. On the other hand , LOW 
corresponds to the matter dominated regime of oscillations in the Earth and LOW 
predicts a well defined structure of peaks in the zenith distribution |34| . The non- 
observation of such peaks in the zenith angle-recoil energy spectrum data decreases 
the likelihood of the upper part of the LOW solution as compared to the LMA 
solution. If one were to use instead the day-night spectra with only night average 
data, this feature would be missed. 

The main difference in the allowed regions shown in the three panels of figure [l] 
is that strategy (b) allows a slightly larger region for the LOW solution. For the 
marginally-allowed (or disallowed) LOW region, the required value for the 8 B neu- 
trino flux can be significantly different from the standard solar model or the SNO NC 
value (the two are virtually indistinguishable). Thus the inclusion of the standard 

3 The strategies (a), (b), and (c) described here correspond, respectively, to the strategies (c), 
(a) , and (b) discussed in detail in ref. Q . 
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solar model uncertainty for the 8 B neutrino flux increases the error used in computing 
the x 2 f° r this strategy, which has the effect of enlarging the allowed region. 

In constructing figure [I], we assumed that only active neutrinos exist. We derive 
therefore the allowed regions in \ 2 using only two free parameters: Am 2 and tan 2 6 4 . 
We use the standard least-square analysis approximation for the definition of the 
allowed regions with a given confidence level. As shown in ref. [^J the allowed regions 
obtained in this way are very similar to those obtained by a Bayesian analysis. 

2.2 Allowed and disfavored solutions 

Table |] gives for our standard analysis strategy (cf. figure [l|a) the best-fit values 
for Am 2 and tan 2 9 for all the neutrino oscillation solutions that were discussed in 
our previous analysis in ref. ETJ. The table also lists the values of Xmin f° r eac h 



,2 



solution. The regions for which the local value of Xmin exceeds the global minimum 
by more than 11.83 are not allowed at 3a CL The number of degrees of freedom in 
this analysis is 46: 44 (Super-Kamiokande zenith-angle energy spectrum) + 2 (Ga 
and CI rates) + 3 ( SNO CC rate , SNO NC rate and A DiV (SNO CC) ) -3 parameters 
( Am 2 , 9, and / B ). 

Within the MSW regime, only the LMA and LOW solutions are allowed at 3a 
with the currently available data. The difference in A% 2 between the global best-fit 
point (in the LMA allowed region) and the best-fit LOW point is A% 2 = 8.8, which 
implies that the LOW solution is allowed only at the 98.8% CL (2.5a). 

The vacuum solution is currently allowed at the 96% CL (2.1a). This solution 
is also allowed in our Before analysis, which we present later in figure || at a CL 
better than 95% CL (and also in the most recent Super-Kamiokande analysis 0). 
But it was not found by the SNO collaboration at the 3a level. In the footnote that 
appears in section |6], we provide a possible explanation for the absence of allowed 
vacuum solutions in the SNO analysis. 

For oscillation solutions for which the survival probability does not depend 
strongly upon energy, such as the LMA and LOW solutions, the effects of includ- 
ing [0, [| the energy distortion in the determination of the rates is very small. We 
find, using the procedure described in section jOj that the central values of the CC 
and NC fitted rates for the best fit points in LMA (LOW) in table p] are shifted by 
+0.5% and —1.5% (0% and +3%), respectively, with respect to the values obtained 
under the hypothesis of no energy distortion. This results in an increase of the A% 2 
between LMA and LOW of ~ 0.2 (i.e., by about 2%). For solutions with stronger 
energy dependences such as SMA (VAC), the effects are larger and lead to shifts in 
the central values of the CC and NC rates of -1.5% and +15.5% (+10% and -11%). 

4 The allowed regions for a given CL are defined in this paper as the set of points satisfying the 
condition 

X 2 (Am 2 ,8) ~xL n < A% 2 (CL, 2 d.o.f.), 
with Ax 2 (CL, 2 d.o.f.) = 4.61, 5.99, 9.21, and 11.83 for CL = 90%, 95%, 99% and 99.73% (3d) 
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Solution 


Am 2 




tan 2 (fl) 


/B,best 


V 2 . 
Amin 


g.o.f. 


LMA 


5.0 


X 


10" 


-5 


4.2 x 10 _1 


1.07 


45.5 


49% 


LOW 


7.9 


X 


10" 


-8 


6.1 x 10" 1 


0.91 


54.3 


19% 


VAC 


4.6 


X 


io- 


10 


1.8 x 10° 


0.77 


52.0 


25% 


SMA 


5.0 


X 


10" 


-6 


1.5 x 10~ 3 


0.89 


62.7 


5.1% 


Just So 2 


5.8 


X 


io- 


12 


1.0 x 10° 


0.46 


86.3 


-0% 


Sterile VAC 


4.6 


X 


io- 


10 


2.3 x 10° 


0.81 


81.6 


-0% 


Sterile Just So 2 


5.8 


X 


io- 


12 


1.0 x 10° 


0.46 


87.1 


-0% 


Sterile SMA 


3.7 


X 


io- 


-6 


4.7 x 10~ 4 


0.55 


89.3 


-0% 



Table 1: Best-fit global oscillation parameters with all solar neutrino data. 

The table gives for the the best-fit values for Am 2 , tan 2 9, Xmin> anc ^ S-°-f- f° r an the 
oscillation solutions among active solar neutrinos that have been previously discussed (see, 
e.g., ref. ||21|). The quantity /b measures the 8 B solar neutrino flux in units of the predicted 
BP00 neutrino flux, see eq. (|2.5| ). The oscillation solutions are obtained by varying the 8 B 
flux as free parameter in a consistent way: simultaneously in the rates and in the night and 
day spectrum fits. The differences of the squared masses are given in eV 2 . The number of 
degrees of freedom is 46 [44 (zenith spectrum) + 4 (rates) + 1 (Adn(CC)) —3 (parameters: 
Am 2 , 6, and /b)]- The goodness-of-fit given in the last column is calculated relative to the 
minimum for each solution. (Solutions that have Xmin — 45.2 + 11.8 = 57.0 are not allowed 
at the 3<T CL) 



We have made a number of checks on the stability of our conclusions regarding 
the global solutions. The A% 2 between LMA and LOW is rather robust under small 
changes in the method of error treatment and in the fitting procedure. However, 
the CL at which the VAC and the QVO solutions (the region between the LOW 
solutions and the VAC solutions, see the insightful discussions by Friedland |35[] and 
Lisi et al. p6|| ) are allowed may fluctuate from just below to just above the 3a limit, 
depending upon details of the analysis. For example, all VAC solutions are disfavored 
at 3a if one ignores the anti-correlation between the statistical errors of the NC and 
the CC rates. 

The best-fitting pure sterile solution is VAC, which is excluded at 5.4<r CL (for 
3d.o.f.). Before April 20, 2002, the best-fitting pure sterile solution was SMA sterile, 
which was acceptable at 3.6cr. 

Oscillations into an admixture of active and sterile neutrinos are still possible as 
long as the assumed total 8 B neutrino flux is increased appropriately ||, 18|. This 



uncertainty can be reduced by combining SNO solar neutrino data with results from 
the terrestrial KamLAND reactor experiment ||. 

Our value for A% 2 = Xlow — Xlma ma Y De lower than most other groups who do 
similar calculations because we include the effect of the correlations and the energy 
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dependence of the neutrino absorption cross sections for the gallium and chlorine 
solar neutrino experiments (see the Appendix of ref. ||). The cross section effects 
correspond to A% 2 = —1.3 for the current global allowed solution, strategy (a). 
Taking this effect into account, our values for A% 2 between the LMA and LOW 



solutions seems to be in general agreement with most other authors [18|, [19|, 20 . 

2.3 Allowed ranges of mass, mixing angle, and 8 B neutrino flux 

The upper limit on on the allowed value of Am 2 is important for neutrino oscillation 
experiments, as stressed in ref. ||22||. In units of eV 2 , we find for the LMA solution 



the following 3a limits on Am 2 

2.3 x KT 5 < Am 2 < 3.7 x 10" 4 . (2.1) 

The LMA solar neutrino region does not reach the upper bound for Am 2 imposed 
by the CHOOZ reactor data WB, i.e., Am 2 < 8 x 10~ 4 eV 2 . Prior to April 20, 2002, 



the solar neutrino data alone were not sufficient to exclude LMA masses above the 
CHOOZ bound. 

For the LOW solution only the following small mass range is allowed, 

3.5 x 10~ 8 < Am 2 < 1.2 x 10~ 7 . (2.2) 



Many authors (see, e.g., ref. []37| and references quoted therein) have discussed 
the possibility of bi-maximal neutrino oscillations, which in the present context im- 
plies tan 2 9 = 1. Figure [I] shows that precise bi-maximal mixing is disfavored for 
both the LMA and LOW solutions. Quantitatively, we find that there there are no 
solutions with tan 2 # = 1 at the 3.3<r CL for the LMA solution, at the 3.2<r CL for 
the LOW solution, and at the 2.8a CL for the VAC solutions. These results refer to 
our standard analysis strategy, corresponding to panel (a) of figure [I]. 

Of course, approximate bi-maximal mixing is now heavily favored. Atmospheric 
neutrinos oscillate with a large mixing angle [Q and all the currently allowed solar 
oscillations correspond to large mixing angles (see figure |IJ). 

How close are the solar neutrino mixing angles to 7r/4? At three sigma, we find 
the following allowed range for the LMA mixing angle 

0.24 < tan 2 9 < 0.89, (2.3) 

and for the LOW solution 

0.43 < tan 2 9 < 0.86. (2.4) 

Let /b be the 8 B neutrino flux inferred from global fits to all the available solar 
neutrino data. Moreover, let /b be measured in units of the best-estimate predicted 
BP00 neutrino flux (5.05 x 10 6 cm" 2 s~ 1 ), 

f - ^ (2 5) 
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The best-fit values for /b for each of the oscillation solutions are listed in the fourth 
column of table [l]. 

The value of /b found by the SNO collaboration from their neutral current 
measurement via a simultaneous solution for all reactions in the SNO detector is [|I] 

/ B = 1.01[1 ± 0.12] (la). (2.6) 

For the global solution shown in figure Qa, the la [3a] allowed range of /b in the LMA 
solution region is 

f B = 1.07 ±0.08 (la) [/ B = 1.07^1] (3a) (LMA). (2.7) 
The 3a range of /b in the LOW solution region is 

/ B = 0.91±g;^ (3a, LOW). (2.8) 



For both eq. ( |2.7| ) and eq. ( |2.8|) , the range of /b was calculated by marginalizing over 
the full space of oscillation parameters (Am 2 , tan 2 28) using the global minimum 
value for \ 2 which lies in the LMA allowed region. 

The allowed range of the 8 B active solar neutrino flux derived from the global 
oscillation solution and given in eq. ( |2.7| ) is slightly more restrictive than was found 
by the SNO collaboration using just their NC measurement (cf. eq. |2.6| ). 

2.4 The Predicted Energy and Day-Night Dependence 

Figure ^| shows the predicted energy dependence and the day-night asymmetry for 
the survival probability P[y e — > v e ) of the allowed MSW oscillation solutions, the 
LMA and LOW solutions, as well as the best-fit vacuum solution. 

We plot in figure |2] the best-fit survival probabilities for an electron neutrino 
that is created in the Sun to remain an electron neutrino upon arrival at a terrestrial 
detector. The right-hand panels present the survival probability for energies between 
and 15 MeV. The left-hand panels present blow-ups of the behavior of the solutions 
at energies less than 1 MeV. 

The energy dependence of both MSW solutions is predicted to be very modest 
above 5 MeV, in agreement with the fact that no statistically significant distortion 
of the recoil energy spectrum has yet been observed in the Super-Kamiokande ex- 
periment [14]. For the SNO CC measurements, table [I| shows that the expected 



distortions of the first and second moments of the electron recoil energy spectrum 
are predicted to be too small to be measurable with statistical significance. 

At energies below 5 MeV, both the LMA and the LOW solutions are predicted 
to exhibit a significant energy dependence. 

At the energies at which water Cherenkov experiments have been possible, the 
survival probability P(v e —>■ v e ) ~ 1/3, i.e., approximately one over the number 
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Figure 2: Survival Probabilities. The figure presents the yearly-averaged best-fit 
survival probabilities determined with strategy (a) for an electron type neutrino that is 
created in the center of the Sun and arrives at a detector on Earth. The full line refers 
to the average survival probabilities computed taking into account regeneration in the 
Earth and the dotted line refers to calculations for the daytime that do not include earth 
regeneration. The dashed line include regeneration at night. The regeneration effects 
are computed for the location of the SNO detector for the right hand panels and for the 
location of the Gran Sasso Underground Laboratory for the left hand panels. There are only 
slight differences between the computed regeneration effects for detectors located at the 
positions of SNO, Super-Kamiokande, and Gran Sasso (see ref. |3J§). The LOW solutions 
in the right-hand panel are averaged over a small energy band, 0.1 MeV, to suppress rapid 
oscillations caused by a sensitive dependence upon the Earth-Sun distance. The vacuum 
solutions are averaged over an energy band of ±0.05 MeV. 

of known neutrinos. According to MSW theory, this is an accident. We see from 
figure || that a survival probability of order 0.5 is predicted for energies less than 1 
MeV. 

The day-night asymmetry is potentially detectable for the LMA solution only at 
energies above 5 MeV; the predicted asymmetry increases with energy (see figure 
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The situation is just the opposite for the LOW solution. The day-night asymmetry 
is large below 5 MeV and relatively small above 5 MeV. 

The vacuum solution has rapid oscillations in energy above 2 MeV (see figure |2]), 
but these oscillations would be very difficult to observe because they tend to average 
out over a typical experimental energy bin. At lower energies, the vacuum oscillations 
will also appear to be smooth, again because of the difficult of resolving in energy 
the oscillations. 

2.5 Analysis details for aficionados: the two-step procedure 

The relatively precise values of the CC and NC rates given in the recent SNO paper [|IJ 
were extracted by a fit to the observed recoil energy spectrum using the assumption of 
an undistorted spectrum, i.e. for a survival probability that is constant in energy 5 . 
These rates cannot be used directly to test an oscillation hypothesis which would 
cause significant distortions of the CC and ES recoil energy spectra. Since they 
were aware of the inconsistency in doing so, the SNO collaboration did not use in 
their oscillation analysis the rates extracted assuming undistorted spectra. Instead, 
they correctly performed a direct fit to their summed spectrum, including the three 
contributions from NC, CC and ES (as well as the background) computed for each 
point in oscillation parameter space. 

We include the spectral distortions in a somewhat different way than was done 
by the SNO collaboration. We perform a two step analysis, which has advantages 
that we discuss at the end of this section. We use the data generously provided by 
the SNO collaboration [||]. 

• First, for a given point in neutrino oscillation parameter space, we compute for 
each of the 17 SNO energy bins the oscillation probabilities that are appropriate 
for evaluating the CC and ES recoil energy spectra. We also compute the 
number of events that would be expected in each bin for an undistorted energy 
spectrum (survival probabilities equal to unity everywhere). We multiply the 
number of events in a given energy bin for the undistorted spectrum by the 
oscillation probabilities and then add the CC and ES energy spectrum to the 
NC energy spectrum (which is undistorted). We fit the observed SNO energy 
spectrum with an arbitrary linear combination of the three energy spectra, 
CC, ES, and NC, plus the background. The coefficients of the three spectra 
are evaluated by minimizing the y 2 fit to the observed SNO spectrum. In 
order to test easily whether the individual results are physically plausible, we 
normalize the coefficients so that they are unity for an undistorted spectrum. 
We have verified that the inferred NC, CC, and ES rates are always within 
the expected physical range. We never encounter unphysical situations: the 

5 They also quote the NC rate obtained from a fit with arbitrary distortion, which results into a 
much less precise determination. 
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extracted rates for all three processes, NC, CC, and ES, are always positive 
and the CC rate never exceeds the NC rate in any part of the parameter space, 
even in those regions with large spectral distortions 6 . In this way, we obtain as 
a function of the oscillation parameters, the best-fit CC and NC measured rates 
(and their corresponding statistical errors, which are strongly anti-correlated) . 

Second, with these NC and CC rates, we compute the \ 2 f° r each particular 
point in oscillation space, combining these two measurements with the avail- 
able results of all other solar neutrino experiments. We use the three different 
analysis strategies discussed in ref. || and in section ^TT] of this paper. We use 
for the SNO day-night asymmetry the result quoted for an undistorted spec- 
trum , since the small effect of the distortion is expected to nearly cancel out 
in the asymmetry. In combining the NC and CC data, we take account of the 
correlation between their errors (both statistical and systematic). We make 
the approximation that the systematic uncertainties have the same percentage 
values [1], 0] as they have for the undistorted rates. This approximation is also 
used, among others, by the SNO collaboration in their recent analysis 0. In 
their recent work, Fogli et.al JiOj discuss the possible effect of this approxima- 



tion and find that the approximation is accurate near the local minima but 
may induce some inaccuracy in the allowed ranges near the boundaries of the 
3a limits. The results of our full-SNO analysis, presented in section. H, confirm 
the conclusion of Fogli et al. We find essentially the same effects as reported 
by Fogli et al. when we no longer make the assumption that the uncertainties 
have the same percentage values as are computed for a undistorted spectrum. 

The constraints on the relations between the ES, CC, and NC rates that exist 
for active oscillations are implemented in step two after the CC and NC ex- 
tracted rates are included in the global fit together with the ES spectra from 
Super-Kamiokande. We do not include the SNO ES rate in our global fit, since 
the Super-Kamiokande ES rate is much more precise, but we do verify, that 
the extracted SNO ES rate is always compatible with the Super-Kamiokande 
ES rate. 

What are the advantages of the two-step procedure? The most obvious advantage 
is speed. It is faster to solve for the allowed neutrino oscillation parameters in \ 2 
space using just the three SNO data points representing the CC and NC fluxes 
and the day-night difference than it is to solve for the allowed parameters using 
all 34 points in the SNO day-night spectrum. Moreover, the use of two distinct 



6 In section 2.2, we noted that the normalized coefficients differ from unity by typically 1% or 2% 
for the LMA and LOW solutions, i. e., the spectra distortion is not important for these cases. For 
the SMA and VAC solutions, which have stronger energy dependences for the survival probabilities, 
the shifts can be ~ 10%. 
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methods, the two-step procedure described here and the full-SNO procedure utilized 
in section [|, permits a test of the robustness of the conclusions. In the two-step 
procedure, the SNO data are represented by only 3 data points in the global x 2 
(compared to two rate parameters for the radiochemical experiments, chlorine and 
gallium). Computation of the global x 2 f° r the full-SNO procedure requires the use of 
34 data points from SNO. Thus for the two-step procedure, SNO is represented in the 
global x 2 by a comparable number of data points as the radiochemical experiments 
while in the full-SNO procedure the SNO experiment contributes more than 10 times 
as many data points as the radiochemical experiments. Both methods, the two-step 
procedure and the full-SNO procedure, give (see section |5.2|) similar allowed regions 
and very similar predicted values for solar neutrino observables that will be measured 
in the future. 

Finally, we note that the two-step method has the advantage of transparency. By 
looking at the computer which contains the coefficients of the NC, CC, and ES energy 
spectra we can immediately see the effect of any choice of oscillation parameters on 
the different event rates. This visual inspection also allows us to check quickly for 
possible unphysical or implausible solutions. 

Given that the spectral energy distortions are strongly constrained by the fact 
that Super-Kamiokande 0, fj] does not observe a significant distortion, it was a 
priori unlikely that the two methods of taking account of distortions-the two-step 
procedure and the full-SNO procedure-would lead to significantly different results 
when analyzing actual solar neutrino experimental data. This expectation is verified 
quantitatively in section |5\2. 



3. Global "Before and After" 

What is the impact of the recent SNO measurements |l], |^] on the globally allowed 
regions of neutrino oscillation parameters? 

Figure |3| compares the allowed regions that are found with the solar neutrino 
data available prior to the presentation of the recent SNO data (i.e., prior to April 
20, 2002) with the allowed regions found including the recent SNO data |], 0. We 
like to refer to figure |3] as our "Before and After" figure. 

The left panel of figure |3| was computed including the improvements that we 
described in section [I] regarding the neutrino cross section errors and correlations 
and the improved average gallium event rate. Hence figure |||a does not correspond 
to any previously published global oscillation solution, although it could have been 
computed prior to April 20, 2002. 

The recent SNO measurements have greatly shrunk the allowed region for the 
LOW solution and have significantly reduced the allowed region for the LMA solution, 
as can be seen by comparing the two panels of figure In particular, maximal mixing 
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Figure 3: Global "Before and After." The left panel shows the allowed regions 
for neutrino oscillations computed using solar neutrino experimental data available prior 
to April 20 2002. The right panel shows the allowed region computed with the same 
procedure but including the SNO NC and improved CC data [jl| and the SNO day-night 
asymmetry Q. In the "Before" panel, the LOW solution is allowed at 97.4% and in the 
"After" panel LOW is allowed at the 98.8%. 

is now not allowed at 3a and the region does not reach the CHOOZ reactor bound, 
Am 2 < 8 x 1(T 4 eV 2 |l|. 

The LOW solution is now allowed at the 98.8% CL Before the recent SNO 
measurements, the LOW solution was allowed at 97.4%. 

4. Predictions for BOREXINO, KamLAND, SNO and a generic 
p — p detector 

We summarize in this section some of the most important predictions that follow 
from the global neutrino oscillation solutions discussed in section [I] and section ^. To 
test the robustness of the predictions given in this section, we have calculated the 
expected values and allowed ranges by two different methods: 1) using the precise 
values for the CC and NC given in ref. [|lj assuming no spectral energy distortions and 
2) taking account of the potential effects of spectral energy distortions on the SNO 



measurements as described in section 2.5. The differences in the calculated values 



and ranges are negligible in all cases. We present here, for consistency with figure [I] 
and our previous discussion, the values calculated taking account of the potential 
spectral distortions in the SNO data. 
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We summarize in section |4.1| the predictions for the MSW solutions. We answer 
the following question in section [4. 2| : How can we distinguish between vacuum and 
MSW solutions? 

4.1 Predictions for MSW solutions 

Table |2| presents the best-fit predictions and the expected ranges for none important 
future solar neutrino observables. The results summarized in table [2] correspond to 
the right-hand side (the "After" panel) of figure [3[ The notation used here is the 
same as in ref. |4j]. 

How much of a difference have the recent SNO measurements made in the ex- 
pected values for future solar neutrino observables? The reader can answer this 
question by comparing table | with table |] . Here, table || presents the best-fit pre- 
dictions and allowed ranges that correspond to using solar neutrino data available 
prior to April 20, 20002. In particular, the values given in table |3] were obtained 
using the global solution shown in the left-hand side, the "Before" panel, of figure [| 

The principal differences between the results summarized in table and table |3] 
are a consequence of the smaller allowed region of the LOW solution shown in the 
"After" panel of figure The changes that result from the modest reduction of the 
LMA allowed region are generally not very significant. 

The predicted day-night asymmetries, An-d, between the nighttime and the 
daytime event rates for SNO and for BOREXINO ( 7 Be v — e scattering detector) 
are given in the first and third rows, respectively, of table § and table [| The results 
presented correspond to an average over one year. The definition of An-d is 

[Night -Day] 
An " D " 2 [Night + Day] ' 

We have used 1 km steps (6371 total points) to compute the values for the day- 
night effect in the SNO CC measurement. Compared to the results obtained using a 
cruder grid of 50 km steps, the best-fit values for the A^_d for the LMA are shifted 
up by ~ 10%. The LOW predictions for A^-d are not significantly affected. The 
main reason for that shift in the values calculated for the LMA is the more accurate 
parameterization of the PREM density EI| of the earth in the external shells. 



The minimum predicted value for the LOW day-night asymmetry, An_e> ( 7 Be), 
in BOREXINO is now a whopping-big 10% at 3a. Prior to April 20, 2002, a 0% 
value of An-d ( 7 Be) was allowed even for the LOW solution. There are no significant 
differences in the Before and After predictions for the SNO CC value of An_d. 

The first and second moments of the SNO CC spectrum for the case of no 
oscillations are (T ) = 7.74MeV and (<7o) = 1-87 MeV. The range of predicted shifts 
with respect to these no-oscillation values is always smaller than ~ 2% within both 
the LMA and LOW 3a regions and is not significantly different from the range found 
in the "Before" analysis. Larger shifts for the moment second moment (7%) are 
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Observable b.f. ±la LMA ±3<r LOW ±3<r 
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Table 2: "After" Predictions. This table presents for future solar neutrino observables 
the best-fit predictions and la and 3a ranges that were obtained by using analysis strategy 
(a) and all solar neutrino data currently available. The best-fit values and uncertainties 
given here correspond to the allowed regions in the right hand panel of figure ||; i.e., this 
table was constructed by including the data made available by the SNO collaboration on 
April 20, 2002 [@, §. The day-night asymmetries for the SNO and the 7 Be (BOREXINO) 
experiments are denoted by An~d and are defined by eq. Q4.1| ). The reduced 7 Be, p — p, 
and 8 B event rates are defined by eq. ( |4.2| ) and eq. (fO|); the reduced KamLAND CC rate is 
defined by eq. ( |4.4j ) . The first moment of the recoil electron energy distribution is denoted 
by 5T for SNO and <5^7 v i s ibie for KamLAND; the second moments are denoted by 5a. The 
threshold of the recoil electron kinetic energy used in computing the SNO observables for 
this table is 5 MeV. For the BOREXINO experiment, we consider electron recoil energies 
between between 0.25 MeV and 0.8 MeV (see ref. |fL6|| ). We present the results for the 
KamLAND reactor observables for two thresholds, Eth = 1-22 and 2.72 MeV. 
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Table 3: "Before" Predictions. This table presents for future solar neutrino observables 
the best-fit predictions and la and 3<r ranges that were obtained by using analysis strategy 
(a) and solar neutrino data available before April 20, 2002. The best-fit and uncertainties 
given here correspond to the allowed regions in the left hand panel of figure ||. The format 
and procedures used in constructing this table are the same as used in constructing table || 
except that here we have not included the recent SNO data [jl], ^]. 



possible within the allowed VAC oscillation region. The non-statistical uncertainties 
in measuring the first and second moments in SNO have been estimated, prior to the 
operation of the experiment, in ref. [12] and are, respectively, about 1% and 2%. 

Based upon the results given in table we predict that SNO will not measure 
a statistically significant (> 3a) distortion to the recoil energy spectrum for the CC 
reaction. This prediction constitutes an important consistency test of the oscillation 
analysis and the understanding of systematic effects in the detector. 
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The prediction for the reduced 7 Be v — e scattering rate, 

7 _ Observed v — e scattering rate , . 

^ ^ = BPOO predicted rate ' ^ ' ' 

is remarkably precise and remarkably stable. The current prediction for what BOREX- 
INO will measure if the LMA solution is valid is [ 7 Be] = 0.64l^g. A somewhat 
smaller value is predicted if the LOW solution is correct. 

We have also evaluated the predictions for a generic p — p neutrino-electron 
scattering detector. The reduced rate for this detector is defined, analogously to the 
reduced rate for 7 Be detectors, by the relation 

r _ Observed v — e scattering rate 

^ ~~ P * = BPOO predicted rate ' ^ ' ' 

We present the predicted rate for two plausible kinetic energy thresholds, 100 keV and 
50 keV. The predicted rate is precise and robust, which supports previous suggestions 
that a measurement of the p — p neutrino scattering rate can be used to determine 
accurately the dominant solar neutrino mixing angle. 

We also include in table |2| the predictions for the expected v — e scattering rate, 
[R( 8 B)] (defined similarly to [ 7 Be] and [p — p}), from 8 B neutrinos above two different 
electron recoil kimetic-energy energy thresholds, T t ^ = 3.5 MeV and T t h = 5 MeV 7 . 
The rates expected if the LMA or LOW solution is correct are essentially identical to 
the neutrino-electron scattering rates measured at higher thresholds by the SNO and 
Super-Kamiokande collaborations, which are 0.47 ±0.13 and 0.46 ±0.04, respectively 
(relative to the value predicted with the BPOO standard solar model flux). Thus a 
measurement of the 8 B v — e scattering rate in BOREXINO, or with a low threshold 
in Super-Kamiokande, is not expected to yield a significant deviation from the rate 
measured at higher energies. 

The predicted value of the reduced CC event rate in the KamLAND reactor 
experiment, 

r^^n/r^ t aatt^n Observed v + p absorption rate , . 

CC KamLAND = -f- ^ , (4.4) 

JNo oscillation rate 

is not significantly affected by the recent SNO results. The best-fit prediction shifts 
slightly (to a lower value) but the shift is well-within the la currently allowed range. 

For KamLAND, it is convenient to represent the distortion of the visible energy 
spectrum by the fractional deviation from the undistorted spectrum of the first two 



moments of the energy spectrum. We follow the notation and analysis of refs. 0, 42 



7 If we normalize the detector exposure to have 55 events a day in the energy window 0.2 MeV < 
T < 0.8 MeV (G. Bellini, private communication, 5/2002) and assume that the detector efficiency 
is energy independent (but include of course the energy dependence of the neutrino cross section 
and the detector resolution of BOREXINO), we estimate for the BPOO v e flux a total of 174 events 
a year above 3.5 MeV of which 109 would be above 5 MeV. 
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43]. The predicted fractional distortion of the first two moments is not significantly 
affected by the recent SNO measurements. Since table || and table || give only 
the fractional changes of the moments relative to the moments for the undistorted 
spectrum, we must also specify the values calculated for a spectrum unaffected by 
new physics. In the absence of oscillations, one expects J4j: (-Evis)o = 3.97 MeV and 
(ct)o = 1.26 MeV for threshold = 1-22 MeV ((E vis ) = 4.33 MeV and (a) = 1.06 
MeV for £ threshold = 1.72 MeV). 

4.2 Predictions for vacuum solutions 

The clearest evidence against the vacuum oscillations would be the observation of a 
rate depletion or a spectral distortion in the KamLAND reactor experiment. The 
Am 2 for vacuum oscillations is too small to lead to an observable effect with Kam- 
LAND. 

Within the 3a allowed VAC regions, the distortion of the SNO spectrum cor- 
responds to a shift in the first (second) moment of the recoil energy distribution of 
at most -2% (+6%) and no significant seasonal variation at SNO is expected. The 
predicted Be rate for vacuum oscillations at BOREXINO is in the same range as 
the predictions for LMA and LOW solutions. The most striking signal for vacuum 
oscillations would be the observation of a large seasonal variation in the BOREXINO 



experiment, with a clear pattern of the monthly dependence of the observed rate [[44 
There should also be a day-night effect at BOREXINO associated with this seasonal 
variation; the day-night asymmetry should be at most ±8% (the size and sign of this 
asymmetry is sensitive to the exact value of Am 2 considered to be within the allowed 
VAC islands) due to the dependence of the survival probability upon the earth-sun 
distance. 



5. Global analysis including the full SNO day-night energy 
spectrum 

In order to help determine the robustness of the oscillation solutions found in sec- 
tion H, we have recalculated the global solutions using the full SNO day-night energy 



spectrum instead of the two step treatment of SNO data described in section g]5 
This full-SNO procedure has been used, by among others, the SNO collaboration [I 



We begin by presenting in section |57J] the global solutions obtained using the full 
SNO analysis procedure and then compare in section [572] the allowed ranges for Am 2 , 



tan 9, and /b found by the two step procedure and the inclusion of the full SNO 



day-night energy spectrum. We compare in section 5.3 the predictions for future 



measurements with the BOREXINO, KamLAND, and SNO detectors and a generic 



future p — p solar neutrino detector. We describe in section |5.4| details of the analysis 



procedure, details that are probably only of interest to aficionados. 
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Figure 4: Global neutrino oscillation solutions with full SNO day-night spec- 
trum and analysis strategy (a). To assess the robustness of the global oscillation 
solutions, compare this figure with figure |]a or ||b. The input data used in constructing 
figure H include the SNO day-night spectrum H, the Chlorine |]l^] and Gallium (aver- 
aged) H [jT|, m [l3|] event rates, as well as the zenith angle-recoil energy spectrum data 
presented by Super-Kamiokande 0] and the neutrino fluxes and uncertainties predicted by 
the standard solar model Q except for the 8 B flux (which is treated as a free parameter). 
The rates from the GALLEX/GNO and SAGE experiments have been averaged to provide 
a unique data point (72.4 ± 4.7 SNU). The CL contours shown in the figure are 90%, 95%, 
99%, and 99.73% (3cr). The global best-fit points are marked by a star. 



5.1 Global solutions 

Figure [| presents the result of our "full-SNO" global analysis (strategy a) that was 
obtained using the full SNO day-night energy spectrum in combination with the 
Super-Kamiokande zenith-angle energy spectrum and the Ga and CI event rates. 
The number of degrees of freedom in this analysis is 77: 44 (Super-Kamiokande 
zenith-angle energy spectrum) + 2 (Ga and CI rates) + 34 (SNO day-night energy 
spectrum) —3 parameters ( Am 2 , 9, and /b)- 

To simplify the comparison with our results obtained with the two-step analysis 
of the SNO data (see figure and figure ||b), we continue to use the average gallium 
rate 72.4 ± 4.7 SNU and do not include the recent, preliminary GNO data reported 
at Neutrino 2002 (which would result into a slight lowering of the rate to 70.8 ± 4.4 
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Solution 


Am 2 




tan 2 (9) 


/B,best 


2 

Xmin 


g.o.f. 


LMA 


5.8 


X 


10 


-5 


4.5 x 10" 1 


1.02 


75.4 


53% 


LOW 


7.9 


X 


io- 


-8 


6.3 x 10" 1 


0.93 


85.0 


25% 


VAC 


6.5 


X 


10" 


-10 


1.5 x 10° 


0.76 


85.5 


23% 


SMA 


4.7 


X 


10 


-6 


1.7 x 10~ 3 


0.92 


101.0 


3.4% 


Just So 2 


5.8 


X 


10" 


-12 


1.0 x 10° 


0.45 


115.1 


0.3% 


Sterile VAC 


4.5 


X 


10" 


-10 


2.5 x 10° 


0.78 


101.4 


3.3% 


Sterile Just So 2 


5.8 


X 


io- 


-12 


1.0 x 10° 


0.45 


111.8 


0.6% 


Sterile SMA 


3.6 


X 


io- 


-6 


4.2 x 10~ 4 


0.55 


115.1 


0.3% 



Table 4: Best-fit global oscillation parameters obtained using the full SNO day- 
night energy spectrum. The values given here, which were computed using the full SNO 
day-night energy spectrum, should be compared with the values given in table [I], which were 
computed using the two-step treatment of SNO data that is described in section |2.5| . Table |] 
gives the best-fit values for Am 2 , tan 2 #, Xmin> an< ^ g-°-f- f° r an the oscillation solutions 
among active solar neutrinos that have been previously discussed (see, e.g., ref. pT|). The 
quantity /b measures the 8 B solar neutrino flux in units of the predicted BP00 neutrino 
flux, see eq. fl2.5|) . The oscillation solutions are obtained by varying the 8 B flux as free 
parameter in a consistent way: simultaneously in the rates and in the fits to the night and 
day energy spectra. The differences of the squared masses are given in eV 2 . The number of 
degrees of freedom is is 77 [44 (Super-Kamiokande zenith-angle energy spectrum) + 2 (Ga 
and CI event rates) + 34 (SNO day-night energy spectrum) —3 (parameters: Am 2 , 9, and 
/b)]- The goodness-of-fit given in the last column is calculated relative to the minimum 
for each solution. (Solutions that have Xmin — 75.4 + 11.8 = 87.2 are not allowed at the 3a 
CL) 

SNU). We have verified that the use of this lower rate would not affect significantly 
any of the conclusions of the present analysis. The x 2 value for LMA would be 
lowered by +0.3, while the difference between LMA and LOW (and VAC) solutions 
is decreased by 0.3. 

5.2 Comparison of two-step and full-SNO analysis procedures 

Table f| gives, for our full-SNO global analysis (cf. figure [|), the best-fit values for 
Am 2 and tan 2 9 for all the neutrino oscillation solutions. The table also lists the 
values of Xmin f° r eacxl solution. The regions for which the local value of Xmin exceeds 
the global minimum by more than 11.83 are not allowed at 3a CL. 

Comparing the right hand panel of figure |3| with figure ^ and the corresponding 
minima in table [l] with those in table |4] shows that the main effect of the inclusion 
of the full day-night spectrum information is, as expected, the moderate worsening 
of solutions with larger spectrum distortions. SMA becomes even further disfavored 
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(acceptable only at the 4.7a CL) while vacuum solutions are now only marginally 
allowed at 3a (acceptable at the 2.7a CL). Furthermore, the best-fit vacuum solution 
is moved into the second lowest Am 2 island. The confidence level for the LOW 
solution exceeds 99% CL (acceptable at the 2.6a CL). 

How are the ranges of allowed masses and mixing changed? 

In units of eV 2 , we find for the LMA solution the following 3a limits on Am 2 , 

2.4 x KT 5 < Am 2 < 4.7 x 1(T 4 . (5.1) 
For the LOW solution, the allowed mass range is 

3.2 x KT 8 < Am 2 < 1.1 x 1(T 7 . (5.2) 



Comparing these results with the corresponding results given in eq. [O] and eq. \L2 
for the two-step procedure for analyzing SNO data, we see that in both cases the 
full-SNO procedure gives a slightly larger 3a range for Am 2 . 
The allowed range for the LMA mixing angle is 

0.27 < tan 2 9 < 0.92 (5.3) 

and for the LOW solution 

0.47 < tan 2 9 < 0.89 (5.4) 

Comparing the limits given above with the values found earlier (see eq. [E3] and 
eq. |2~4p using the two-step procedure, we conclude that the full-SNO procedure gives 
a slightly shifted 3a range for tan 2 9. 

For the global solution shown in figure |], the lcr[3cr] allowed range of f B is 

f B = 1.02 ± 0.08 (la) [f B = 1.02 ±0.24] (3a) (LMA). (5.5) 

The range of f B was calculated, as just as it was for eq .(|2.7|) and eq. (|2.8|), by 
marginalizing over the full space of oscillation parameters (Am 2 , tan 2 29). Hence 
we use Ax 2 conditions for 1 dof, relative to the global minimum value for \ 2 that 
lies in the LMA allowed region. Since for the full-SNO analysis procedure A% 2 = 
Xlow ~ Xlma > 9, there are no allowed solutions for / B at the 3<r CL (1 dof) within 
either the LOW or vacuum solution domains. 

The uncertainties shown in eq. |5.5| are essentially the same as the uncertainties in 
/ B that were found using the two-step procedure for analyzing SNO data (see eq $L7\ 
and eq. |2.8| ). We conclude that the uncertainties in all three of the parameters Am 2 , 
tan 2 9, and f B are robust with respect to the choice of analysis procedures for the 
SNO data. 
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5.3 Predictions for BOREXINO, KamLAND, SNO and a generic p-p 
detector 

The predicted event rates for measurable solar neutrino quantities are essentially 
the same if one uses the full-SNO day-day night energy spectrum or if one uses the 
two-step procedure. 

Table || gives the predicted ranges for future solar neutrino observables that 
correspond to the allowed regions in figure |]. The values shown in table |^ should be 
compared with the results presented in table |], which was obtained using the two-step 
procedure for SNO data. The comparison shows that there are small shifts in some 
of the central values of the observables due to the slight shifts in the positions of the 
LMA and LOW best fit points. However, the predicted ranges are not significantly 
modified. 



5.4 Analysis details for aficionados: using the full SNO day-night energy 
spectrum 

We describe in this subsection the details of our analysis procedure when we include 
the full SNO day-night energy spectrum. We perform a direct fit to the total SNO 
day-night energy spectrum, including contributions from NC, CC, and ES events. 

In order to avoid redundancy with the discussion in section p~5| , we also provide in 
this section some detailed information about analysis procedures and error estimates 
that is necessary for comparing our results with the results obtained by other analysis 
groups. This additional information is relevant to both the two-step and the full-SNO 
analysis procedures. 

We use the detailed information provided in ref. || by the SNO collaboration. 
For each point in oscillation parameter space, we compute (as with the two-step 
procedure) the expected event rates for the chlorine and gallium experiments, and 
the Super-Kamiokande zenith-angle energy spectrum. For the SNO day-night energy 
spectrum, we compute the NC, CC, and ES contributions in each energy bin. We 
construct a theoretical SNO spectrum by summing the NC, CC, and ES contributions 
together with the background rates given in ref. ||. With these model predictions, 
we make a global fit including as the 80 data points, the 2 radio-chemical rates, the 
44 data points of the Super-Kamiokande zenith-angle spectrum, and the 34 points 
of the SNO day-night spectrum. 

We include the errors and their correlations for all of these observables. We 
construct an 80 x 80 covariance error matrix that includes the effect of correlations 
between the different errors as off diagonal elements (an alternative "pull" method 



has been recently proposed in ref. [40]). The main sources of correlations are due 
to the BP00 neutrino fluxes, the cross section uncertainties for the gallium and 
chlorine rates, and the cross sections and experimental systematic errors for Super- 
Kamiokande and SNO energy spectra. As in our two-step analysis, the flux and cross 
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Observable 


b.f. ±lcr 


LMA ±3(7 


LOW ±3cr 


A N _ D (SNO CC) (%) 


4.2 + nq 


4.2+I-2 

— 4.Z 


2.7 + o'o 


XT WMn cc\ fo/\ 
01 (^olNU KjKj) y/o) 


-0-10_o.5o 


-0.10_ 164 


0.46_ 041 


5a (SNO CC) (%) 


U ' ZO -0.50 


9 25 +0 - 39 


0-38l°;g 


[R ( 7 Be)] v — e scattering 


0.64 ±0.03 


64+ - 09 


0.59 ±0.04 


A N _ D ( 7 Be) (%) 


o o+°-° 

u - u -o.o 


+ai 


23i? 3 


[T? f^~R^l 7/ £> Cr*Q ^"i"OT*lTl o* 

IV I J-J J\ Is O oCdb LCI J-lJ-g, 








(T th = 3.5 MeV) 


0.46 ±0.03 


0.46 ± 0.04 


0.45 ±0.03 


(T th = 5 MeV) 


0.45 ±0.03 


0.45 ± 0.04 


0.45 ±0.03 


fCCl (KamLAND 
(E th = 2.72 MeV) 
(E th = 1-22 MeV) 


0-58±g:i? 
0.571J" 


c i8 +0 - 14 
S7 +a12 

u - d ' -0.30 


— 


5E ■ u, (KamLAND"! (%) 

v -'--'visible ^J-YdiiiJ-j-riiN u j J 

(E t h = 2.72 MeV) 
(E th = 1-22 MeV) 


1 —2 

p;+10 
~~ °-5 


-7+i 4 

1 —4 

p;+13 
~~ °-8 


— 


(5(7 fKamLAND) f%) 
(E th = 2.72 MeV) 
(E th = 1-22 MeV) 


-Hlf 


-11+f 

-14±g 2 




\p-p] v — e scattering 








(T th = 100 keV) 
(T th = 50 keV) 


0.699™ 
0.695l°;°| 


0.699±85SS 
u.uyj_ 4 5 


0.6861^1 
0.6811^1 



Table 5: "After" Predictions obtained using the full SNO day-night energy 
spectrum. This table presents for future solar neutrino observables the best-fit predictions 
and lex and 3<r ranges that were obtained by using the analysis including the SNO day- 
night spectrum and all other solar neutrino data currently available. The best-fit values and 
uncertainties given here correspond to the allowed regions in the figure ||. The predictions 
given in this table should be compared with table 0, which was obtained by treating the 



SNO data by the two-step procedure described in section 2.5. 



section errors are included following the general approach described in ref. |3(J with 
the refinements described in the appendix of ref. H. 

The experimental systematic errors can be either correlated or uncorrected in 
energy. We include these errors as follows 8 . For Super-Kamiokande, we take from 



8 We describe here our error treatment for both Super-Kamiokande and SNO energy spectra. 
For Super-Kamiokande, the treatment is the same for the two-step analysis and for the full-SNO 
procedure. 
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table 1 of ref. |7j the systematic error, which we include as uncorrelated among the 
Super-Kamiokande energy bins but fully correlated among the zenith-angle bins. For 
SNO, energy independent but correlated errors include the systematic errors from 
the uncertainty in vertex reconstruction (+3% for CC and ES and +1.45% for NC), 
the background errors from neutron capture, and low-energy correlated systematic 
uncertainties. Each of these errors are fully correlated among the 34 SNO bins. 
Furthermore, the fractional error is assumed to be constant, i.e. independent (in %) 
of the oscillation point and affecting equally 8 B and hep neutrinos. 

For both Super-Kamiokande and SNO there also energy-dependent and energy- 
correlated errors from (i) the 8 B spectral energy shape, (ii) the absolute energy 
scale, and (iii) the energy resolution. The fractional value of those errors is, in 
general, affected by distortions of the spectrum due to oscillations and therefore 
they should be evaluated for each point in the oscillation parameter space. The 
Super-Kamiokande and the SNO collaborations do not specify the details of what 
they use for these errors. We evaluate the three errors listed above as follows. First, 
we compute for each point in oscillation parameter space the expected number of ES 
events in each energy bin for the Super-Kamiokande spectrum, and the number of 
CC, ES, and NC events in each energy bin for SNO. We then evaluate the fractional 
change in the calculated results by: (i) shifting the normalization of the 8 B spectrum 
up and down within the la error determined in ref. j|5] in order to obtain the 
error associated with the 8 B spectrum shape; (ii) shifting the calibrated energy T' to 
T'(1±S T ) with 5T = 0.64% (1.21%) for Super-Kamiokande (SNO) to obtain the error 
associated with the energy scale; (iii) sifting the width of the resolution function a to 
o-(l ± 8 a ) with 5 a = 0.025 for Super-Kamiokande(5 CT = 0.045 + 0.004 x (T — 4.98) for 
SNO) to obtain the error due to the energy resolution. Whenever a comparison was 
possible, we have verified that the error estimates made as described above are in 
good agreement with the information about the errors that has been provided by the 
Super-Kamiokande and SNO collaborations. Each of these errors is fully correlated 
among the different bins in each experiment. Moreover, the shape errors are also 
correlated between the Super-Kamiokande and SNO experiments. 



6. Discussion and summary 

In the following, we present without parentheses results that were obtained by the 
two-step procedure for analyzing SNO data (and present in parentheses the corre- 
sponding results that were obtained by analyzing the full SNO day-night energy 
spectrum). If the results are identical for both treatments of the SNO data, then 
we omit the value in parenthesis. By comparing the values given with and without 
parentheses, the reader can see the rather small differences that result from different 
choices of how to treat the SNO data. 
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Figure [I] and figure [| together with table |] (or, alternatively, figure || and ta- 
ble (|), tell much of the story. These global analyses, which use all the available 
solar neutrino data, demonstrate that only large mixing angle solutions are currently 
allowed at the 3.7a (4.7a) confidence level. We therefore know at an impressive 
confidence level that there are relatively large mixing angles for oscillations between 
solar neutrinos (and between atmospheric neutrinos), unlike the small mixing angles 
among quarks. However, precise maximal mixing is excluded at 3.2a for MSW so- 
lutions [see eq. |2l| and eq. |2]4| (eq. ^]3| and eq. |5]4]) for the exact allowed regions] 



and at 2.8a (2.9a) for vacuum oscillations. The masses and mixing angles for solar 
neutrinos do not appear to satisfy an obviously simple pattern. 

Among the MSW solutions, the situation has also been clarified. The LMA 
solution is now the only viable solution at a level of 2.5a (2.6 a). The LOW solution 
is excluded at the 98.8% (99.2%) CL. The KamLAND and BOREXINO experiments 
will test strongly this conclusion. The SMA solution is now excluded at more than 
3.7a (4.7 a). Pure sterile oscillations are excluded at 5. 4a (4.7 a). 

The possibility of a strictly energy independent solution of the solar neutrino 
problem is excluded at more than 3.6a. If we consider solutions that have an energy 
variation less than 10% over the energy range of interest, the exclusion is at more than 
2.6a. We note that such a weak energy dependence can probably not be distinguished 
experimentally from strict energy independence (see ref. ||6| for details). 

Vacuum solutions are not favored, but are acceptable at 2.1a (2.7a) 9 . 

For the reader who (like us) often prefers pictures to tables, figure ^ shows the 
predicted energy dependence and day-night difference for the current best-fit LMA, 
LOW, and VAC solutions. 

The global oscillation solutions constrain the active 8 B solar neutrino flux with 
somewhat greater accuracy than the SNO NC measurement. The accuracy of the 



9 We have attempted to determine why the SNO collaboration did not find any 3a allowed 
vacuum solutions. We have found at least one plausible explanation. If we make the imprecise 
approximations that are listed below, then we also do not find any allowed vacuum solutions at 
the 3cr CL. Here are the approximations (with references to descriptions of more accurate error 
treatments) : 1) neglect the uncertainty in the true shape of the 8 B undistorted energy spectrum, 
2) neglect the dependence on the neutrino oscillation parameters of the calculated errors for the SNO 
and Super-Kamiokande spectral energy distributions (i. e., compute the errors for an undistorted 
8 B energy spectrum so that the percentage errors for each energy bin are independent of Am 2 
and tan 2 8 ), and 3) neglect the correlations and energy dependence of the errors |, fF§ for the 
chlorine and gallium neutrino cross sections. Making these approximations, we find that the A% 2 
between the best-fit LMA solution and the vacuum solution with the lowest \ 2 is 14.6 for the 
two-step SNO procedure and is 12.3 for the analysis including the full SNO day-night energy 
spectrum. Both of these values exceed the A\ 2 of 11.83 that corresponds to a 3tr CL. We also 
tested the effect of neglecting the energy shape |Q of the 7 Be line. This approximation turned out 
to be important for predicting the values of future observables for BOREXINO and KamLAND, 
but docs not significantly affect the computed difference between the global best-fit values of x 2 
between LMA and vacuum solutions. 
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Source 


CI 


Ga 


CI 


Ga 


CI 


Ga 




(SNU) 


(SNU) 


(SNU) 


(SNU) 


(SNU) 


(SNU) 




LMA 


LMA 


LOW 


LOW 


VAC 


VAC 


DD 





40.4 





38.2 





40.3 


pep 


0.12 


1.51 


0.10 


1.25 


0.15 


1.82 


n on 
11CJJ 


n m 


n D9 




n rn 

u.uo 


n D9 




7 Be 


0.62 


18.6 


0.53 


16.0 


0.46 


14.0 


8 B 


2.05 


4.35 


2.26 


4.72 


2.34 


5.00 


13 N 


0.04 


1.79 


0.04 


1.56 


0.05 


1.83 


15Q 


0.15 


2.83 


0.15 


2.44 


0.18 


3.01 


17p 


0.00 


0.03 


0.00 


0.03 


0.00 


0.04 



Total 3.03 69.6 3.11 64.2 3.21 66.0 

Table 6: Neutrino oscillation predictions for the chlorine and gallium radio- 
chemical experiments. The predictions are based upon the global analysis strategy (a), 
and use the neutrino fluxes for the BP00 standard solar model (except for the 8 B flux for 
which the factor /b is included) and the neutrino absorption cross sections H, The 



rates are presented for the best-fit oscillation parameters of the allowed solutions listed in 
table |l|. The total rates should be compared with the standard solar model values ||, which 
are, 7.6^}'^ (chlorine) and 128^ (gallium), and the measured values, which are 2.56 ± 0.23 
(chlorine @) and 72.4 ± 4.7(gallium g 0, ||. ) 



SNO NC measurement is ±12% [|T|. We find from a global oscillation solution of all 
the data, that at the la CL, / B = 1.07 ± 0.08 (/ B = 1.02 ± 0.08), where f B is the 
total flux of active 8 B solar neutrinos in units of the flux predicted by the standard 
solar model ||. The allowed ranges of /b for the LMA and LOW solutions are given 
in eq. (|2.7|) and eq. fl2~8| ) (eq. [531) at la and 3a CL. 

We have calculated the predicted best-fit values and allowed ranges for ten im- 
portant solar neutrino observables for the BOREXINO, KamLAND, and SNO ex- 
periments and for a generic detector of v — e scattering by p — p neutrinos. The 
observables include total rates, day-night effects, and spectral distortions. The pre- 
dictions are summarized in the "Before" and "After" tables, table |2| and table |^, and 
in section 12. Similar results are presented for the full-SNO procedure in table (H) 
and in section |5.3[ The predictions are rather robust and are not significantly af- 
fected by the choice of procedure (two-step or full-SNO) used to analyze the SNO 
data. Only the predictions for the disfavored LOW solution are significantly affected 
by the recent SNO measurements |], ||. We have also calculated the predicted val- 
ues and ranges given in table and table |3] both with and without taking account 
of potential spectral energy distortions in interpreting the SNO data (see section 
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and ref. 0). Both methods give essentially identical results. 

For completeness, table |6] gives the individual contributions of different solar 
neutrino sources to the predicted event rates for the chlorine and gallium experiments. 
The results presented in table |6] were obtained using the best-fit oscillation solutions 
of table |I[ The values given in the table were calculated with the two-step treatment 
of the SNO data. The calculated chlorine and gallium rates are essentially the same 
when the best-fit solutions are used that were obtained including the full SNO day- 
night energy spectrum. In that case, the predicted rates for the LMA solution are 
increased by about 1% and the predicted values for the LOW and vacuum solutions 
are decreased by between 2% and 7%. The rates shown in table || are also very 
similar to the rates that were expected |4]] prior to the recent SNO measurements. 

The most important next generation solar neutrino experiment will detect the 
fundamental p-p neutrinos, which constitute 91% of the total neutrino flux predicted 
by the BPOO solar model. Unfortunately, there are no approved p-p solar neutrino 
experiments at the present time, although there are a number of very promising 
proposals under development. Among the CC detectors, only LENS |jD] is in a ad- 
vanced stage of research development and even for LENS important details regarding 
the detector are not yet known. There are a number of potential v — e scattering 
experiments |)T| that could observe p-p neutrinos. Therefore, for completeness we 
have added the predictions for a generic v — e scattering experiment using p-p neu- 
trinos. Table || and table |3] (see also table |^) give the calculated values of the p — p 
v — e scattering rate for two plausible recoil electron kinetic energy thresholds, 50 
keV and 100 keV. 

We are grateful to the SNO collaboration for the thrill of analyzing solar neutrino 
data that include a neutral current measurement. We have profited from several 
discussions with M. Chen, F. Duncan, and A. McDonald regarding the analysis 
methods used by the SNO collaboration and from a constructive referee's report. 
JNB acknowledges support from NSF grant No. PHY0070928. MCG-G is supported 
by the European Union Marie-Curie fellowship HPMF-CT-2000-00516. This work 
was also supported by the Spanish DGICYT under grants PB98-0693 and FPA2001- 
3031. 
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